Utility Map

#acs_vars <- load_variables(2022, "acs5")

chas_df <- load_chas() %>%
    mutate(State_Name = case_when(
    state == 17 ~ "Illinois",
    state == 25 ~ "Massachusetts",  
    state == 27 ~ "Minnesota",  
    state == 34 ~ "New Jersey",
    state == 35 ~ "New Mexico",  
    state == 36 ~ "New York",
    state == 06 ~ "California",
    TRUE ~ "Other"
  )) %>%
  filter(State_Name == "Minnesota" )

mn_chas <- chas_df %>%
  filter(State_Name == "Minneosta")

get_tracts <- get_acs(
  geography="tract", 
  state="MN",
  variables=c("total_pop" = "B01001_001", 
              "total_hhs" = "B11001_001", 
              "ami_total" = "B06012_001",
              "ami_one" = "B06012_002",
              "ami_two" = "B06012_003"),
  year=2019, 
  geometry=TRUE) %>%
  group_by(geoid=GEOID) %>%
  summarize(total_pop = estimate[variable=="total_pop"],
            total_hhs = estimate[variable=="total_hhs"],
            pop_below_150 = estimate[variable=="ami_one"] + estimate[variable=="ami_two"],
            ami_total = estimate[variable=="ami_total"])

## Getting data from the 2015-2019 5-year ACS

## Downloading feature geometry from the Census website.  To cache shapefiles for use in future sessions, set `options(tigris_use_cache = TRUE)`.

add_geo <- merge(chas_df, get_tracts) 

add_geo <- st_as_sf(add_geo)